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Abstract 

Coronaviruses are single-stranded positive-sense RNA viruses associated with important avian diseases. Their relatively high rates of 
mutation and recombination frequencies allow them to adapt to new hosts and ecological niches. Although Brazil has 18% of global 
avian species diversity, studies regarding the presence of avian viral diseases in wild birds in South America are scarce. In this study, 
we performed a retrospective analysis of the presence of CoVs in 746 wild birds. Oropharyngeal and cloacal swabs were obtained and 
placed together in vials containing VTM transport medium collected in different regions of Brazil between 2006 and 2013. Screening 
for viral nucleic acid was performed using conventional RT-PCR and pancoronavims nested PCR. Positive samples were character¬ 
ized by partial sequencing of the RNA-dependent RNA polymerase (RdRp) gene, and ensuing phylogenetic analysis was performed 
to investigate the association between vims epidemiology and bird migration routes. Coronavims RNA were detected and sequenced 
from six samples, in which three were related to gammacoronavimses group and the other three to deltacoronavims group. Our study 
documents the presence of Co Vs related to avian gamma- and deltacoronavimses circulating in both urban- and poultry-farm regions 
of Brazil, implicating wild birds as potential carriers of Co Vs which may represent a risk to poultry farms and public health in Brazil. 
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Introduction 

Coronavimses are RNA single-stranded positive polarity vims 
associated with important disease in avian poultry [1]. 


Infectious bronchitis vims (IBV), known as coronavims of 
chickens, causes important economic losses to the poultry in¬ 
dustry. Studies about the presence of Co Vs in wild birds, es¬ 
pecially IBV, in South America are scarce [2]. Coronavimses 
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(CoVs) possess the largest genomes among RNA viruses and 
have relatively high mutation rates and recombination frequen¬ 
cies when compared with other single-strand RNA viruses, 
which may allow them to adapt to new hosts and ecological 
niches [3-5]. These characteristics are exemplified by the two 
emerged highly pathogenic betacoronaviruses, CoV-SARS 
(severe acute respiratory syndrome) and the CoV-MERS 
(Middle East respiratory syndrome) [6]. 

Even though avian coronavimses with zoonotic potential have 
not been reported yet, many species of wild and domestic birds 
can carry Co Vs responsible for substantial economic loss within 
the poultry industry such as the gammacoronavims infectious 
bronchitis vims (IBV). This vims affects the gastrointestinal, 
respiratory, urinary, and reproductive systems and may lead to 
mortality rates as high as 30% and sudden drops in egg produc¬ 
tion [7]. Moreover, many novel Co Vs genetically distinct from 
IBV have been recently identified in different families of wild 
and domestic birds and mammals but their potential economic 
losses remain unknown. Recently, Duraes-Carvalho et al. (2015), 
described the detection of delta- and betacoronaviruses in captive 
wild birds, from a rehabilitation center, originally captured in Sao 
Paulo and Mato Grosso do Sul state with close evolutionary 
relationships with Co Vs isolated from mammals [8]. 

Wild migratory birds may facilitate the spread of vimses, as 
demonstrated by the surveillance of avian influenza vims [9, 
10] and West Nile vims [1 1-1 3], but the circulation of CoV in 
avian migratory routes in South America, more specifically in 
Brazil, still need to be elucidated. The potential diversity of 
CoV in Brazilian bird fauna as consequence of the bird species 
diversity (18%), together with the large number of endangered 
species and the high intensity of poultry farms in Brazil, indi¬ 
cates that birds from Brazilian biomes are important subjects 
for CoV surveillance [2, 14]. 

In this study, we performed a retrospective analysis, mo¬ 
lecular detection, and phylogenetic analysis of avian 
coronavimses using oropharyngeal and cloacal swab samples 
from wild and domestic birds collected between 2006 and 
2013 in different regions of Brazil to investigate the associa¬ 
tion between vims epidemiology, bird migration routes, and 
the proximity of urban areas and poultry farms. 

Materials and methods 
Sample collection 

In total, 746 samples (Table 1) of birds captured between 2006 
and 2013 in five expeditions to north, northeast, midwest, south, 
and southeast Brazil (Fig. 1) were analyzed. For selection, the 
migratory routes were investigated through a map analysis and 
experience of ornithologists [15]. The samples were collected by 
field team of the Laboratory Clinical and Molecular Virology, 
University of Sao Paulo (USP), Sao Paulo, Brazil, with support 


of the Wildlife Department of Sao Paulo (DEPAVE); University 
of Vale do Rio dos Sinos (UNISINOS); University of 
Pernambuco (UFRPE); National Center for Research and 
Conservation of Wild Birds (CEMAVE); Laboratory of Bird 
Ecology, Institute of Biosciences (UFMT); and Laboratory of 
Virology and Rickettsioses, Veterinary Hospital (UFMT). 
Samples were collected from wild birds after capture with mist 
nets, strategically placed between 30 and 250 cm of the ground, 
and by hand nets on an open field in public parks. The field 
works were made according to the permission of the Brazilian 
Institute of Environment and Renewable Natural Resources 
(IBAMA) and National Center for Research and Conservation 
of Wild Birds (CEMAVE) developed and followed protocols 
approved by the Ethics Committee on Animal Experimentation 
of the Institute of Biomedical Sciences, University of Sao Paulo. 
Sampling from oropharyngeal and cloacal swabs was obtained 
and placed together in vials containing VTM transport medium 
[16]. All birds were then released after sampling. 

RNA extraction and reverse transcription 

Viral RNA was extracted using the MagMax TM-96 RNA 
Isolation Kit (Ambion, Austin, TX, USA) according to the 
manufacturer’s instructions. The extracted RNA was eluted 
in 80 pL of RNase-free water and stored at - 80 °C until 
processed. The cDNA was transcribed using the High- 
Capacity cDNA Archive Kit (Applied Biosystems, Foster 
City, CA, USA) following the manufacturer’s instructions. 
The RT reaction was performed by adding, to each reaction, 
10 pL of total RNA, 2 pL of RT Buffer 10X, 2 pL of Random 
Primers (50 ng/pL) (Thermo Fisher Scientific, Waltham, MA, 
USA), 4.2 pL of DEPC water, and 50 U of MultiScribe™ 
Reverse Transcriptase (Thermo Fisher Scientific). The RT step 
conditions were 25 °C for 10 min, 37 °C for 120 min, and 
85 °C for 5 min, and cooling to 4 °C. 

Coronavirus detection 

Polymerase chain reaction (PCR) followed by nested PCR 
were performed according to the protocol for the detection 
of avian Co Vs developed by Chu et al. (2011) [17], which 
targets the RNA-dependent RNA polymerase (RdRp) gene 
in the ORF1 area of the viral genome. The resulting approx¬ 
imately 600-bp PCR fragment and 440-bp nested PCR frag¬ 
ment were separated by agarose gel [1.5%] electrophoresis 
and viewed under UV light after staining with ethidium bro¬ 
mide (25 pg/100 mL). 

Purification and sequencing 

Samples were purified using the ExoSap-IT enzyme (GE 
Healthcare, Chalfont St. Giles, UK) and GeneJET™ Gel 
Extraction Kit (Thermo Fisher Scientific) according to the 
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Table 1 Surveillance of Coronavirus in resident and wild birds in Brazil by molecular tests. Bold represents positive samples 


Site 

Date 

Scientific name 

Sample* 

No. pos./No. birds 

Ilha de Marajo, PA 

Jul. 2006 

Cairina mochata 

C and O 

0/42 



Galius gallus 

C and O 

0/36 



Meliagres gallopavo 

C and O 

0/15 

Ajuruteua, PA 

Nov.2006 

Actitis macularius 

C and O 

0/5 



Calidris pusilla 

C and O 

0/2 



Charadrius semipalmatus 

C and O 

0/1 



Dendrocygna autumnalis 

C and O 

0/18 



Arenaria interpres 

C and O 

0/9 



Dendrocygna viduata 

C and O 

0/56 

Ilha de Canelas, PA 

Nov. 2006 

Actitis macularius 

C and O 

0/10 



Calidris canutus 

C and O 

0/1 



Calidris pusilla 

C and O 

0/3 



Calidris minutilla 

C and O 

0/1 



Sterna hirundo 

C and O 

0/4 



Larus sp. 

C and O 

0/1 


Nov. 2008 

Actitis macularius 

C and O 

0/31 



Arenaria interpres 

C and O 

0/22 



Calidris alba 

C and O 

0/2 



Limnodromus griseus 

C and O 

0/9 



Tringa melanoleuca 

C and O 

0/1 



Charadrius semipalmatus 

C and O 

0/1 



Gelochelidon nilotica 

C and O 

0/1 

Bahia de Sao Marco, MA 

Apr. 2007 

Dendrocygna viduata 

C and O 

0/27 



Netta erythrophthalma 

C and O 

0/2 



Dendrocygna autumnalis 

C and O 

0/48 



Gallinula chloropus 

C and O 

0/2 



Amazonetta brasiliensis 

C and O 

0/14 



Anas bahamensis 

C and O 

0/13 



Meliagres gallopavo 

C and O 

0/1 

Lagoa do Peixe National Park, RS 

Nov. 2009 

Rynchops niger 

C and O 

1/63 



Calidris fuscicolis 

C and O 

1/63 



Haematopus palliatus 

C and O 

0/6 



Himantopus melanurus 

C and O 

0/2 



Pluvialis dominica 

C and O 

0/7 



Arenaria interpres 

C and O 

0/4 



Calidris canutus 

C and O 

0/3 



Calidris alba 

C and O 

1/6 



Calidris pusilla 

C and O 

0/1 



Nycticriphes semicollaris 

C and O 

0/1 



Fulica leucoptera 

C and O 

0/1 



Pujfinus gravis 

C and O 

0/1 



Pujjinus puffinus 

C and O 

0/3 



Tringa flavipes 

C and O 

0/2 



Sterna hirundo 

C and O 

0/15 



Phalacrocorax brasilianus 

C and O 

0/2 



Vanellus chilensis 

C and O 

0/1 



Elaenia mesoleuca 

C and O 

0/1 



Pluvialis squatarola 

C and O 

0/1 



Tryngites subruficollis 

C and O 

0/7 



Springer 





Braz J Microbiol 


Table 1 (continued) 

Site 

Date 

Scientific name 

Sample* 

No. pos./No. birds 



Calidris melanotus 

C and O 

0/1 



Sterna vitatta 

C and O 

0/1 

Sao Jose do Egito Caatinga, PE 

Oct. 2010 

Thaminophilus capstratus 

C and O 

0/8 



Hylophilus amaurocephalus 

C and O 

0/10 



Sarkesphorus cristapus 

C and O 

0/8 



Hemitricus margaritaseiventer 

C and O 

0/4 



Coryphospingus pileatus 

C and O 

0/10 



Anopetia gounellei 

C and O 

0/1 



Casiornis fuscus 

C and O 

0/1 



Gyolophylox hellmayri 

C and O 

0/3 



Tolmyias flaviventris 

C and O 

0/2 



Heliomaster squamosus 

C and O 

0/2 



Cychlars gujanensis 

C and O 

0/1 



Cyanoloxia brissoni 

C and O 

0/1 



Formicivora melanogastrer 

C and O 

0/8 



Cantorchilus longirostris 

C and O 

0/1 



Chlorostilbon lucidus 

C and O 

0/1 



Coereba flaveola 

C and O 

0/1 



Polioptila plumbea 

C and O 

0/1 



Synalax frontalis 

C and O 

0/1 

Pantanal Pocone, MT 

Nov. 2012 

Ramphocelus carbo 

C and O 

0/10 



Hem itriccus m argarita ceiven ter 

C and O 

0/6 



Cnemotriccus fuscatus 

C and O 

0/11 



Paroaria capitata 

C and O 

0/2 



Picumnus albosquamatus 

C and O 

0/3 



Gallus gallus 

C and O 

0/9 



Furnarius rufus 

C and O 

0/8 



Myiarchus swainsonii 

C and O 

0/2 

Ibirapuera Park, SP 

Jun. 2013 

Anser cygnoides 

C and O 

3/81 

Total 


63 species 


6/746 


* Sample type: C = cloacal swabs and O = oral swabs 


manufacturer’s instructions. The Sanger sequencing reaction 
was prepared using Big Dye Terminator v3.1 Cycle 
Sequencing Kit (Life Technologies, Foster City, CA, USA), 
followed by purification using BigDye X Terminator® 
Purification Kit (Applied Biosystems), and analysis in an au¬ 
tomated ABI PRISM 3130XL DNA Sequencer (Applied 
Biosystems). 

Genetic characterization and phylogenetic tree 
of the RdRp gene 

Generated chromatograms containing nucleotide sequences 
were visualized using Bioedit Sequence Alignment Editor 
(v.7.1.5.0). Nucleotide sequences were aligned using CLC 
sequencer WorkBench version 7.6 (QIAGEN Aarhus, 
Denmark), and the nucleotide percentage similarities were 


calculated using the MegAlign program of the DNASTAR 
package (DNASTAR, Inc., Madison, WI, USA). The similar¬ 
ities with existing sequences were determined using the Basic 
Local Alignment Search Tool (BLAST). Phylogenetic analy¬ 
sis was performed using the neighbor-joining method (NJ) 
implemented in the MEGA 7.0.26 [18]. The nucleotide sub¬ 
stitution model used for phylogenetic reconstruction was the 
p-distance model. The bootstrap values >55% were obtained 
in the analysis of 10,000 replicates and are presented at the 
branching points. 


Fig. 1 Location of sampling sites with positive results for avian CoVs ► 
and their relationship with the main bird migratory routes and 
congregation sites in North and South America. Congregation sites 
include reproduction, wintering, and tagging areas of boreal migrant 
birds (modified from http://www.icmbio.gov.br/cemave/) 


<£) Springer 






Braz J Microbiol 



9*ntirva 


TWfridtlFMgo 


Hwdiion Bay 


tates of Ame 


Unit 

L‘. J'** 




Legend 

Migratory Route 
Reproduction Site 
Wintering Site 
Staging Site 
^ Sampling Site With Positive Results 

Sampling Site Without Positive Results 


wmwm 

F _ 

A Delaware Bay 


• irm 


K Njlha da Marajd* PA 
PA ' ISha de Canelas* PA 

Maranhao 


Pantanal-MT 


SJ Panta 

•ll 


> f 7 2 


doAvlAo-Pf 


Maogoe 

Se<o* BA 


Sio Paulo- SP 


P Lagoa do Wi# 
National Parle* RS 


San Antonio Oette 


4^ Springer 



Braz J Microbiol 


Table 2 Sampling location, total 
samples, number of positive 
samples obtained from birds 
collected across Brazil from 2006 
to 2013 


Location 

No. of samples 

Positive samples 

Year 

Ilha de Marajo/PA 

93 

— 

2006 

Ajuruteua/PA 

91 

— 

2006 

Ilha de Canelas/PA 

67 

— 

2006-2008 

Bahia de Sao Marco/M A 

107 

— 

2007 

Lagoa do Peixe National Park/RS 

192 

3 

2009 

Sao Jose do Egito/PE 

64 

— 

2010 

Pantanal Pocone/MT 

51 

— 

2012 

Ibirapuera Park/SP 

81 

3 

2013 

Total 

746 

6 

2006-2013 


PA Para, SP Sao Paulo, MT Mato Grosso, BA Bahia, RS Rio Grande do Sul 


Results 

The CoV screening protocol, based in a nested PCR assay 
with degenerate consensus primers targeting the conserved 
region of RNA-dependent RNA polymerase (RdRp) gene, 
was selected for its greater sensitivity and CoV genera ampli¬ 
tude detection. Of the 746 bird samples collected from 2006 to 
2013 in different regions of Brazil and analyzed by conven¬ 
tional RT-PCR, six samples (0.8% of total) from two different 
regions were confirmed by DNA sequencing (Table 2). Our 
study describes the detection of samples related to avian 
gamma- and deltacoronaviruses in birds collected in southeast 
and southern region of Brazil, more specifically in Sao Paulo 
and Rio Grande do Sul state (Table 3). The presence of CoV 
RNA was detected in three Chinese geese (Anser cygnoides ) 
sampled in Ibirapuera Park, downtown Sao Paulo, southeast 
Brazil. One sample has shown similarities (93%) with se¬ 
quences from viruses in the gammacoronavirus group, which 
are related to migratory birds only, whereas the other two 
samples (84%) belong to the deltacoronavirus group. Even 
though these sequences are closely related to those of wild 


bird Co Vs, this group also includes ferret, pig, and leopard 
Co Vs, demonstrating the potential for spillover events be¬ 
tween coronaviruses from bird and mammal hosts. In the same 
way, in southern Brazil, three samples positive for Co Vs were 
obtained from migratory birds captured in wintering areas at 
Lagoa do Peixe National Park, state of Rio Grande do Sul. 
Phylogenetic analysis showed that these samples (PNLP 100, 
PNLP 115, and PNLP 159) belonged to two CoV genera. 
Samples related to gammacoronavimses were obtained from 
one Calidris alba and one Calidris fuscicollis individual, re¬ 
spectively, related a restricted to sandpiper species (Fig. 2). 
Samples related to deltacoronavirus were obtained from a 
black skimmer (Rynchops niger) carrying a strain that clus¬ 
tered, in phylogenetic analysis, closer to the clade of viruses 
found in Ibirapuera Park (Fig. 3). 

Discussion 

The nested PCR assay with degenerate consensus primers was 
used to many studies. This assay has been used for the 


Table 3 Sample identification, bird species, place of capture, and coronavirus group of avian coronavirus sequences obtained from bird samples 
collected across Brazil from 2006 to 2013 


Sample 

Species 

Location 

Group 

Year 

Accession 

number 

(GenBank reference) 

Similarity (%) 

DPV 5 

Anser cygnoides 

Ibirapuera Park/SP 

Deltacoronavirus 

2013 

KU321643 

Wigeon coronavirus HKU20 
strain HKU20-9243 (2012) 

84 

DPV 10 

Anser cygnoides 

Ibirapuera Park/SP 

Deltacoronavirus 

2013 

KU321644 

Wigeon coronavirus HKU20 
strain HKU20-9243 (2012) 

84 

DPV 16 

Anser cygnoides 

Ibirapuera Park/SP 

Gammacoronavirus 

2013 

KU321645 

Avian coronavirus isolate 

7 M ZPorphyrula alleni (2011) 

93 

PNLP 100 

Calidris alba 

Lagoa do Peixe 
National Park/RS 

Gammacoronavirus 

2009 

KU321640 

Rock Sandpiper Coronavims 
(2010) 

99 

PNLP115 

Rynchops niger 

Lagoa do Peixe 
National Park/RS 

Deltacoronavirus 

2009 

KU321641 

Chroicocephalus 

ridibundus/Finland (2013) 

81 

PNLP 159 

Calidris fuscicollis 

Lagoa do Peixe 
National Park/RS 

Gammacoronavirus 

2009 

KU321642 

Western Sandpiper 

Coronavirus/KR-28 (2010) 

99 


SP Sao Paulo, RS Rio Grande do Sul 
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Lagoa do Peixe National Park 
| Ibirapuera Park 

I Positive Control- IBV 


C/5 

> 

03 

C 

o 

o 

o 

03 

0 

Q 


0.06 


-cE 



GQ427173.1 /Turkey/coronavirus/T CoV/VA-74/03 
J EU095850.1 /Turkey/coronavirus/MGI 0 

- KJ690953.1 /Mallard/coronavirus/ACoV/Mallard/PL-MW345/2009 
KJ135013.1/Chicken/infectious_bronchitis_virus/IBVUkr27-11 
Positive Control-Avian Infectious bronchitis virus strain Ma5 
KJ425488.1/Chicken/infectious_bronchitis_virus/ck/CH/LHB/110825 

KJ690955.1/Bean_Goose/coronavirus/ACoV/Bean_goose/PL-MW435/2009 

KM093893.1/Porphyrula_alleni/coronavirus/CV11/2011 
KM093897.1 /Bulbucus _ibis/coronavirus/6cvPF 12/2011 
KM093895.1 /Gallinago_macrodactyla/coronavirus/4cvPF7/2011 
KM093899.1 /Bulbucus _ibis/coronavirus/9cvPF 17/2011 
KM093879.1 /Charadrius_pecuarius/coronavirus/S22/2011 
KM093885.1 /Gallinula_chloropus/coronavirus/4P/2011 
KM093883.1/Charadrius_pecuarius/coronavirus/2P/2011 
KM093882.1 /Rostratula_benghalensis/coronavirus/1 P/2011 
KM093881.1/Gallinula_chloropus/coronavirus/S23/2011 
KM093880.1 /Anas_hottentota/coronavirus/S21 //2011 
PNLP100/Calidris_alba/coronavirus/Brasil/2009 
95 J PNLP159/Calidris_fuscicollis/coronavirus/Brasil/2009 
GU396675.1 /Western_sandpiper/coronavirus/KR-28 
GU396687.1/Rock_sandpiper/coronavirus/CIR-65824 
L GU396685.1/Rock_sandpiper/coronavirus/CIR-665821 

KM093892.1/Porphyrula_alleni/coronavirus/7M/Porphyrula_alleni/2011 

KM093890.1/Porphyrula_alleni/coronavirus/8CV/Porphyrula_alleni/2011 

KM093901.1/Bulbucus_ibis/coronavirus/6cvPF11/Bubulcusibis/2011 

KM093891.1/Porphyrula_alleni/coronavirus/5P/Porphyrula_alleni/2011 

DPV_16/Anser_cygnoides/coronavirus/Brasil/2013 

58 — GU396690.1/Snow_goose/coronavirus/WIR-159 

I GU396677.1 /Brent_goose/coronavirus/KR-69 

L GU396676.1 /Brent_goose/coronavirus/KR-70 

— JN788879.1/DendrocygnaJavanica/coronavirus/KH08-0852/Dendrocygna 


Fig. 2 Phylogenetic analysis of partial RdRp gene sequence of avian 
gammacoronaviruses. Phylogenetic tree of partial sequences of the 
RdRp-encoding region constructed using the neighbor-joining method. 
In green, CoVs related to gammacoronaviruses sequences detected in 
samples from Lagoa do Peixe National Park, state of Rio Grande do 


Sul, southern Brazil. In red, CoVs related to gammacoronaviruses 
sequences detected in samples from Ibirapuera Park, Sao Paulo, southeast 
Brazil. Bootstrap values > 55% were obtained in the analysis of 10,000 
replicates and are presented at the branching points 


identification and differentiation of CoV in the four genera 
described to date (alpha-, beta-, gamma-, and 
deltacoronavirus) [17, 19, 20]. Similarly to our report, most 
of avian CoV diversity studies detected the same coronavirus 
genera [2, 17, 21, 22], with the exception of a report of 
betacoronavirus detection in distinct birds species from a re¬ 
habilitation center in Brazil [8]. Further studies are needed to 
understand the circulation of betacoronavirus in birds. 

Ibirapuera is the main park in Sao Paulo, receiving on 
average 150,000 visitors in a single weekend. Wild birds 
make up a large portion of the species of animals that in¬ 
habit the park (ca. 142 species). The birds are kept free and 
cohabiting with the human population and other wild species 
(resident or visitors). There are another 21 species of wild 
animals, including bats and skunks, and the park also allows 
pets [23]. Such close contact between so many different 
species may cause interspecific transmission and 


recombination events, leading to the appearance of emerging 
viruses [24]. The three goose samples were collected from 
geese that live in this same environment in contact with wild 
birds that use the park as a stopping point. On the other 
hand, three migratory birds captured in wintering areas at 
Lagoa do Peixe National Park-RS have shown positive re¬ 
sults of CoVs too. Sandpipers are boreal migrants that mi¬ 
grate south from North America in the winter [25]. A CoV 
related to deltacoronavirus group was obtained from a black 
skimmer (Rynchops niger ), thus highlighting the possibility 
that migratory birds may be carrying these viruses across 
migratory routes in Brazil [9]. The Lagoa do Peixe region 
is an important migratory route [25], and it is located close 
(400 miles) from one of the most important region of poultry 
farm production of Brazil, more specifically region of Porto 
Alegre in Rio Grande do Sul State [26]. The proximity be¬ 
tween the migratory routes of wild birds infected by 
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100 


58 


100 


86 


58 


81 


KJ601779.1/Swine/Deltacoronavirus/PDCoV/USA/lllinoisl 36/2014 


KJ481931.1/Swine/Deltacoronavirus/PDCoV/USA/lllinoisl 21/2014 
ioo KP981395.1/Porcine/Deltacoronavirus/USA/IL/2014/026PDV 
KJ769231.1/Swine/Deltacoronavirus/OhioCVM1/2014 
KJ601780.1/Swine/Deltacoronavirus/PDCoV/USA/Ohiol 37/2014 

KJ567050.1/Porcine/Deltacoronavirus/8734/USA-IA/2014 
KM012168.1 /Porcine/coronavirus/HKU 15 
EF584908.1 /Leopard_cat/coronavirus/Guangxi/F230/2006 
KP757890.1/Porcine/Deltacoronavirus/CHN-AH-2004 
EF584911.1/Chinese_ferret_badger/coronavirus/Guangxi/F247/2006 
JQ065042.1 /Porcine/coronavirus/HKU 15 
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KJ584356.1 /Porcine/coronavirus/HKU 15 
KM820765.1/Porcine/Deltacoronavirus/KNU 14-04 
JQ065045.1 /Sparrow/coronavirus/HKU 17 


ioo 


FJ376619.2/Bulbul/coronavirus/HKU 11 -934 
— FJ376620.1/Bulbul/coronavirus/HKU11-796 
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FJ376621,1/Thrush/coronavirus/HKU 12-600 

JQ065044.1 /White-eye/coronavirus/HKU 16 
JX548304.1/Shorebird/coronavirus/85 
100 KX588674.1/Chroicocephalus_ridibundus/Finland/2013 

PNLP115/Rynchops_niger/coronavirus/Brasil 

-JQ065049.1/Common-moorhen/coronavirus/HKU21 


ioo 


FJ376622.1/Munia/coronavirus/HKU 13-3514 
JQ065046.1 /Magpie-robin/coronavirus/HKU 18 

DPV_5/Anser_cygnoides/coronavirus/Brasil 


ioo 
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DPV_10/Anser_cygnoides/coronavirus/Brasil 
JQ065048.1/Wigeon/coronavirus/HKU20 


Fig. 3 Phylogenetic analysis of partial RdRp gene sequence of avian 
deltacoronaviruses. Phylogenetic tree of partial sequences of the RdRp- 
encoding region constructed using the neighbor-joining method. In green, 
CoVs related to deltacoronavirus sequences detected in samples from 
Lagoa do Peixe National Park, state of Rio Grande do Sul, southern 


Brazil. In red, CoVs related to deltacoronavirus sequences detected in 
samples from Ibirapuera Park, Sao Paulo, southeast Brazil. Bootstrap 
values >55% were obtained in the analysis of 10,000 replicates and are 
presented at the branching points 


gammacoronavirus shows the importance of a continuous 
surveillance and potential contact control between poultry 
and migratory birds as measures to avoid any chance of viral 
transmission between these animals. 

The CoVs related to deltacoronaviruses group (DPV 5 and 
DPV 10) and gammacoronavirus group (DPV 16) from 
Ibirapuera Park showed the highest similarity to a CoV found 
by Lau et al. (2012) [27] in Guangdong province (84%), 
China, and by Lima et al. (2015) [28] (93% nt identity) in a 
migratory bird (Porphyrula alleni) from Madagascar in 2011, 
respectively (Table 2). The gammacoronaviruses found at 
Lagoa do Peixe National Park were more similar to those 
found by Muradrasoli et al. (2010) [29] in sandpipers 
0 Calidris spp.) from the Bering Strait (99% nt identity). 
However, the deltacoronavims was more similar to the strain 
isolated from a black-headed gull (Chroicocephalus 
ridibundus) by Hepojoki et al. (2017) in Finland (95% nt 
identity) (Table 3). Interestingly, this data were reported about 


sampling collected during the same years (2010-2013) by the 
same way, using RT-PCR targeting the RdRp gene to detect 
the presence of CoV in wild birds from Finland. The preva¬ 
lence was found around 5.4% (51/939), which corroborates 
with our findings. Hepojoki et al. (2017) have shown positive 
samples to diverse gamma- and deltacoronavims too [22]. 

In South America, there is a lack of information about 
coronavirus in wild birds. These data are interesting, once 
that, present birds infected in central park in Sao Paulo 
city and wild birds through migratory routes. In Hong 
Kong, Chu et al. (2011) detected a high prevalence 
(12.5%) in aquatic wild birds, reporting gamma and 
deltacoronavims, being the first one more prevalent [17]. 
Similar studies screened before in Sweden (2016), with 
the number of sampling similar to our study, and reported 
a high prevalence (18.7%) in wild birds, with emphasis 
for Anseriformes that had the highest prevalence (51%), 
all belonging to the gammacoronavirus group [30]. On the 
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other hand, Jordan et al. (2015) analyzed 303 samples 
from cloacal swabs, using universal CoV primers, and 
found only one duck coronavirus, belonging to the 
gammacoronavirus group too, which suggested a low 
prevalence in wild aquatic birds in USA, corroborating 
with our data [21]. As we know, migratory routes in 
Brazil are linked with flyways coming from North 
America and can confirm our findings [16]. Although 
our data have shown low prevalence when compared with 
studies before in Scandinavian Waterfowl, our data is 
clear when compared with North American data. In addi¬ 
tion, it is the first time that samples related to Gamma and 
deltacoronavirus group in free wild birds were reported in 
Brazil. 

The overall prevalence of CoV found in this study was 
0.8% (6/746 birds), with 3.7% (3/81) from Ibirapuera Park 
and 1.6% (3/192) from migratory birds captured in Lagoa do 
Peixe National Park. Other studies report an overall prevalence 
of CoV in aquatic and migratory birds between 0.33 and 18.5% 
[8, 17, 21, 30, 31]. This data reinforces that CoVs are circulat¬ 
ing across major migratory routes worldwide. 

Our report points to the circulation of coronaviruses 
among migratory and resident birds sampled from distinct 
regions of Brazil. Coronavirus RNA were detected in birds 
species of Anatidae family (Chinese goose) in Ibirapuera 
Park, Sao Paulo, southeast Brazil. Anatidae birds were al¬ 
ready described as hosts of coronavirus all over the globe. 
However, interestingly, we report for the first time the detec¬ 
tion in two sandpiper species (Calidris alba and Calidris 
fuscicollis ) from Lagoa do Peixe National Park, a wintering 
area in southern Brazil. Although our data indicate a low 
prevalence of coronavirus in birds sampled, the genetic sim¬ 
ilarity between coronavirus detected at birds from distinct 
continents may indicate that wild birds are possibly carrying 
coronavirus throughout migratory routes as observed for avi¬ 
an influenza virus. We report the detection of samples related 
to delta- and gammacoronaviruses in different species of 
birds, some sampled in a park located inside the largest city 
in South America, indicating the circulation of coronavirus in 
urban areas (Table 3). In addition, all birds tested in this study 
were in apparently healthy populations. Further studies are 
required to determine the species and full genome similarity 
with other avian CoVs. However, our findings indicate that 
Co Vs may be circulating among one of the major avian mi¬ 
gratory routes in the world. 
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